Magnetically-driven electronic phase separation in the semimetallic ferromagnet EuB 6 
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From measurements of fluctuation spectroscopy and weak nonlinear transport on the semimetal- 
lic ferromagnet EuBe we find direct evidence for magnetically-driven electronic phase separation 
consistent with the picture of percolation of magnetic polarons (MP), which form highly conducting 
magnetically-ordered clusters in a paramagnetic and 'poorly conducting' background. These differ- 
ent parts of the conducting network are probed separately by the noise spectroscopy/nonlinear trans- 
port and the conventional linear resistivity. We suggest a comprehensive and 'universal' scenario for 
the MP percolation, which occurs at a critical magnetization either induced by ferromagnetic order 
at zero field or externally applied magnetic fields in the paramagentic region. 



The observation that the magnetic state of a sys- 
tem critically affects its electronic transport properties 
is at the heart of spintronics research. Fundamentals 
of this effect can be studied in a number of more or 
less complex model systems showing large negative (or 
even colossal) magnetoresistance (MR) behavior, e.g. 
magnetic semiconductors and mixed-valence perovskites. 
These materials often have rich phase diagrams, in which 
many phases exhibit intrinsic, i.e. non-chemical, elec- 
tronic phase separation. Nanoscale phase separation in 
turn is thought to play a critical role in the emergence 
of colossal MR (CMR) in such materials Therefore, 
electronic phase separation has been a subject of inten- 
sive recent theoretical and experimental interest. For the 
present study, we have chosen the low-carrier density fer- 
romagnetic (FM) semimetal EuBg. EuBg has a cubic 
lattice symmetry and is magnetically isotropic (Eu 2+ lo- 
calized 4f spins in a 8 SV/2 Hund's rule ground state), yet 
it shows interesting physics, where the interplay of metal- 
licity and the formation of clustered magnetic phases can 
be studied on the fundamental level in a 'clean' system. 

Despite the simple lattice and magnetic structure, the 
physical properties at low temperatures, in particular 
the mechanism of FM ordering in EuBg and its inter- 
play with CMR behavior is far from being fully under- 
stood. While the material undergoes the paramagnetic 
(PM) to FM transition, it exhibits two anomalous fea- 
tures at T C1 ~ 15.5 K and T C2 ~ 12.5 K in electronic 
transport and specific heat measurements, which initially 
had been interpreted as different kinds of FM ordering 
[H [3J. Applying small magnetic fields drastically sup- 
presses the resistivity at the higher transition T Cl . Differ- 
ent mechanisms have been discussed in order to explain 
such a CMR effect in nonmanganite systems, e.g. the sup- 
pression of critical magnetic fluctuations with externally 
applied magnetic fields 0] or a derealization of carri- 



ers due to the overlap of magnetic polarons (MP) [5H7]. 
MP, which first have been suggested in an experimen- 
tal study of the magnetic semiconductors Eui^Gd^Se 
in 1967 jS], are formed when it is energetically favorable 
for the charge carriers to spin polarize the local moments 
over a finite distance, i.e. the localization length of the 
charge carriers [5]. The size of the polaronic clusters is 
thus determined by the balance of the increase in kinetic 
energy of the charge carriers due to their localization and 
the reduction of exchange energy due to alignment of the 
local moments |10j . From small angle neutron scattering 
experiments, De Teresa et al. have demonstrated the ex- 
istence of such clustered phases above the FM ordering 
temperature Tc of manganite materials 5j . Although the 
formation of MP has been discussed for various different 
magnetic systems, the underlying microscopic nature of 
electrical transport in the MP phases is not yet properly 
understood and is a matter of current debate [H [TTHT3] . 

In EuBg, magnetic phase separation has been sug- 
gested, e.g. from muon spin rotation experiments [14] . 
However, direct evidence for electronic phase separation 
and percolation from transport measurements has been 
lacking so far. Yet, an important hint came from recent 
Hall effect measurements by Zhang et al., who observed 
a distinct change in the slope of the Hall resistivity in the 
PM phase (THUS]. The switching field in the Hall effect 
depends linearly on temperature and extrapolates to the 
paramagnetic Curie temperature of the material. The 
authors interpret the switching field occurring at a sin- 
gle critical magnetization as the point of percolation for 
patches of a more conducting and magnetically-ordered 
phase in a PM background. This picture yields excellent 
scaling in an empirical two-component model |16j . pro- 
viding a measure of the degree of electronic phase sepa- 
ration. 

Motivated by these findings, we performed fluctua- 
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tion spectroscopy and weak nonlinear transport (third- 
harmonic voltage) measurements, which — unlike the 
conventional linear resistivity — are sensitive to the mi- 
crogeometry of the sample, i.e. intrinsic electronic in- 
homogeneities. We find direct evidence for electronic 
phase separation and magnetically-driven percolation, 
both when cooling through the FM transition and when 
applying magnetic fields in the PM temperature regime. 

Single crystals of EuB 6 were grown from Al flux as de- 
scribed in [T7]. Fluctuation spectroscopy measurements 
were carried out using a standard ac method in a four- 
terminal setup [151 [E] ■ The (linear and nonlinear) trans- 
port measurements were carried out in standard four- 
probe geometry using ac lock-in technique at a frequency 
of 17 Hz. For weakly nonlinear transport, current density 
and electrical field are related by j = (<j + b\E\ 2 )E, where 
a is the Ohmic or linear conductivity and b the nonlinear 
conductivity coefficient, with b\E\ 2 < a [Ml HQ- The 
cubic nonlinearity can be accessed in ac transport mea- 
surements by detecting a third-harmonic voltage signal 
V^, see e.g. [20l [22j [23] . In general, noise spectroscopy 
and nonlinear transport are a measure of the fourth mo- 
ment of the current distribution |22| and therefore probe 
the microgeometry of the electronic system, which is not 
accessible by the linear (Ohmic) resistance, being related 
to the second moment of the current distribution |24j . 
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FIG. 1. (Color online) (a) Resistivity of the EuB6 sample 
(#1) used for the noise measurements shown in Fig.[5] Inset: 
Enlarged low-temperature region. Peaks in the temperature 
derivative dp/ 'AT mark the two transitions at T C1 = 15.2 K 
and T C2 = 12.6 K. (b) The coefficient of weak nonlinear 
transport measured by third-harmonic voltage generation, 
k = piu/Puj = Vzw/V w , for sample #2. 



FigureJTJa) shows the bulk resistivity of a representa- 
tive EuBg sample (#1, optimized shape for noise mea- 
surements). Upon cooling, the resistivity decreases from 
room temperature down to T* ~ 35 K, where a broad 
minimum is observed. Below that temperature, the re- 
sistivity increases and goes through a maximum at about 
~ 16 K, before it rapidly decreases. As seen in the inset 
of Fig.jlja), below this sharp drop there is a shoulder 



marking the second, lower transition. In the literature, 
the two transitions are often defined as pronounced peaks 
in the temperature derivative dp/dT observed here at 
T Cl = 15.2 K (just below the resistivity maximum) and 
T C2 = 12.6 K, see inset of Fig.[I|a). Such clearly visi- 
ble transitions are observed only for high-quality sam- 
ples [25]. Interestingly, T Cl and T C2 coincide with the 
onset and peak, respectively, of a weak nonlinear con- 
tribution to the electronic transport, shown in Fig.Qb), 
measured by the ratio of third-harmonic and linear resis- 
tivity P3ui/puj shown here for sample #2. Clearly, in this 
temperature region the electronic system exhibits sub- 
stantial local inhomogeneities. We note that although 
there is a slight sample-to-sample variation in T Cl and 
T C2 (16.1 K and 12. 4K, respectively for sample #2, and 
15 K and 12.3 K for #3), the relation between the kink in 
resistance and the peak in the nonlinear signal is ubiqui- 
tous in the samples we have measured. 
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FIG. 2. (Color online) (a) Temperature dependence of the 
normalized noise PSD S R /R 2 (T, f = 1 Hz). Inset exemplarily 
shows a typical spectrum obtained at T = 13 K, being close 
to Sr oc 1/f (solid line). The solid line in the main panel 
shows a fit parameterizing a diverging behavior with a peak 
at T = 12.4 K and a width AT = 2.4 K. (b) Same data in a 
log- log plot highlighting a step-like increase of the 1//- noise 
level below T* ~ 35 K and a power-law divergence following 
roughly Sr/R 2 oc T~ g (solid line), (c) Scaling Sr/R 2 oc R w 
in the percolation regime T C2 < T ~ T C1 . A fit to the data 
yields w ~ —2.1 (solid line). 



Next, we focus on the measurements of the resistance 
noise power spectral density (PSD), Sn(f,T). A typical 
spectrum, measured at T = 13 K, is shown in the inset 
of Fig.[2ja). In the entire temperature range we observe 
l//-type noise, characterizing the intrinsic equilibrium 
conductance fluctuations of the sample. As required, 
this excess noise is quadratic in the current through the 
sample (not shown). The main panel of Fig.[2]ja) shows 
the normalized resistance noise PSD Sr/R 2 {T) taken at 
1 Hz. The data can be described by a parameterized di- 
vergence with the peak at a temperature T = 12.4 K 
coinciding with T C2 and a width AT = 2.4 K. Further- 
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more, two striking features in the noise measurements are 
(i) a step-like increase below about 35 K coinciding with 
the minimum of the resistivity at T* (seen more clearly 
in Fig.[2j[b)), and (ii) a power-law divergence starting 
above T Cl with a pronounced peak at T C2 , below which 
the noise level rapidly drops again. As discussed in de- 
tail below, we interpret this behavior as a reflection of 
the percolation of the poorly conducting component of 
the sample, with the peak determining the percolation 
threshold at T C2 . Between T C2 and T Cl , the noise scales 
as Sr/R 2 oc T 13 with /3 ~ —6.5 ± 0.6, remarkably sim- 
ilar to the power law observed when approaching the 
metal-insulator transition in bulk samples of disordered 
P-doped Si [26]. 

In general, diverging l//-type resistance (or conduc- 
tance) fluctuations are typical for percolative metal- 
insulator transitions [37J [55] ■ In a classical percolation 
scenario, the resistance and resistance noise, respectively, 
are expected to scale as R, Sr/R 2 oc (p — p c ) _t ' K . Here, 
p is the fraction of unbroken bonds of a random resis- 
tor network, p c the percolation threshold, and t and k 
the critical exponents for the resistivity and noise, re- 
spectively. Since the microscopic details determining p 
are not accessible in most cases, it is convenient to ex- 
press the normalized noise as a function of resistance 
Sr/R 2 oc R w , with w = n/t [29 . The scaling exponents 
t, k and w depend on the type of percolating system and 
its dimensionality and are determined by computer sim- 
ulations, where exact solutions are not known [30]. The 
divergence of the noise close to p c usually is explained 
by the reduced number of effective current paths, which 
results in the suppression of cancellation of uncorrclatcd 
resistance fluctuations along different paths, which are 
abundant far away from p c . This effect is similar to the 
inverse dependence of the noise PSD on the volume of 
macroscopic bodies |31j . In a composite system, p may 
be the portion of the metallic phase in a less conduct- 
ing background, the portion of which is q. In certain 
cases, the percolation threshold p c of that metallic phase 
can be different from the percolation threshold q c of the 
less conducting or insulating portion [2T1[32]. Therefore, 
depending on which contribution dominates, the perco- 
lation threshold probed by noise spectroscopy/nonlinear 
transport on the one hand and conventional linear re- 
sistivity on the other hand may be different. In EuBg, 
we consider the MP as entities of a more conducting 
and magnetically-ordered phase in a paramagnetic and 
'poorly conducting' background [7j [16] , which form links 
at T Cl corresponding to percolation threshold p c and 
therefore a continuous conduction path through the sam- 
ple leading to derealization of holes and hence the drop 
of the resistance. Indeed, Siillow et al. [Jj have ascribed 
T Cl = Tm to a metallization transition via the overlap of 
magnetic polarons. The separation of the charge dereal- 
ization and bulk magnetic ordering at T C2 = Tq [7J then 
implies electronic and magnetic phase separation. The 



noise measurements clearly show a percolation thresh- 
old q c at T C2 , which means that the resistance and its 
fluctuations are sensitive to different parts of the con- 
ducting network. Consequently, we observe a negative 
scaling exponent w ~ —2.1 ± 0.3 above the peak of the 
noise, i.e. for q > q c , where the resistivity increases with 
increasing temperature. Remarkably, such a clear sepa- 
ration is not a unique property of EuBg but also has been 
observed in perovskite manganites [33], where, however, 
a two-component percolation scenario is not discussed. 
Close to the threshold, the continuum 3D 'random void' 
and 'inverted random void' models have scaling expo- 
nents w — 2.1 and 2.4, respectively [34]. For the present 
data, the scaling on the low-temperature side, i.e. below 
the percolation threshold, could not be verified due to 
the limited number of data points. 

As shown in Fig.[l](b), the characteristic temperatures 
T Cl and T C2 are clearly seen also in the weak nonlinear 
transport, as the onset and pronounced peak, respec- 
tively. In a percolative system close to the threshold, 
the local current densities and electric fields can be much 
larger than their average values. This is due to narrow 
paths, so-called 'bottlenecks' and 'hot spots' in the ran- 
dom network, the contribution of which to the resistance 
is proportional to the local electric field squared. Since 
R{E) « R u + AE 2 + ■ ■ ■ with the second term being 
the third-harmonic nonlinear resistance, the magnitude 
of i?3 W indirectly reflects the density and local arrange- 
ment of the MP, see also [23 . Therefore, we consider the 
electronic system as consisting of magnetically ordered 
spheres, which — when they overlap — form more con- 
ducting clusters, characterized by a conductivity tr + , em- 
bedded in less conducting regions with o~ . Both Sr/R 2 
and P3u;/Pai peak at T C2 , the percolation threshold q c for 
the less conducting regions, which is reached when the 
MP clusters merge and start forming a continuum. Most 
interestingly, T C2 (and q c ) coincides with the material's 
bulk Curie temperature as determined from magnetiza- 
tion measurements [7J. This important observation may 
be understood by taking into account that the zero-field 
magnetic susceptibility diverges at Tc, which is most fa- 
vorable for the stabilization of MP. Hence, at this tem- 
perature, the density of MP should be maximum, which 
is the criterion for the percolation threshold we observe. 
In theoretical studies it is found [2U] that weak nonlinear 
transport and 1// noise behavior of a composite system 
of conducting spheres in a less conducting background, 
such as the present one, is more complex than for classi- 
cal bond or site percolation. This is because of the role 
of the local microgeometry, parameterized by the chan- 
nel width between the conducting spheres and the degree 
of their overlap, which in turn is related to substantially 
enhanced local electric fields. These local geometric pa- 
rameters, the degree of disorder of the MP (i.e. deviations 
from a regular array) and the ratios cr + /er~ and S^/S^ 
determine whether the good or the bad conductor dom- 
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inates the overall noise and the nonlinear transport sig- 
nal. In the region of the FM transition, obviously the 
contribution of the background electrons, i.e. the poor 
conductor, dominates. 
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FIG. 3. (Color online) (a) R^ as a function of magnetic field 
H measured at different temperatures in the PM regime. In- 
set shows the same data normalized to the field value and 
magnitude of the maximum in a log-log plot, demonstrating 
that the data collapse on one curve [35 • (b) Phase diagram 
of EuB6 determined from the field sweeps in (a). Circles de- 
note the position of the maximum in R^u (H) showing a linear 
behavior, (c) Schematic of the MP states at different tem- 
peratures in zero magnetic field (after [T3]). Note that the 
schematic is a 2D projection of spherical MP in 3D. 



The MP are suggested to form (to become stabilized) 
below T* ~ 35K [1 [3 [H]. At about the same tem- 
perature, the noise shows a pronounced step-like in- 
crease upon cooling, which can be interpreted as hopping 
of magnetic clusters, the number of which increases as 
the temperature is lowered, and/or scattering by these 
magnetic excitations. The observed peak structure in 
Sr/R 2 {T) may reflect the energetics of such a mecha- 
nism. In the PM region T Cl < T < T* , the MP are ini- 
tially isolated and diluted in the electronic 'background 
sea'. When a magnetic field is applied, the MP increase in 
size [36] until they form links, eventually overlap (where 
at the same time their number may decrease) and fi- 
nally merge. Figure[3ja) shows a pronounced maximum 
in i?3^(i/) (corresponding to q c ) for different T in the 
PM regime, which succeeds the field-induced percolation 
of the MP (corresponding to p c ), similar to the situa- 
tion when cooling in zero magnetic field through the FM 
transition. Indeed, a drop of the sample resistance pre- 
cedes the peak in (not shown). A striking observa- 
tion is a linear temperature dependence of Rf£ x (H), see 
Fig.^b). Thus, like the switching field in the Hall effect 
[16) , the maximum in nonlinear transport occurs at a sin- 
gle critical magnetization. The magnitude of i?3 [i) (T, H) 



is related to the microgeometry of the phase-separated 
electronic system reflecting the density and arrangement 
of the MP. Accordingly, at high temperatures, slightly 
below T* , only a few MP are stable and large fields are 
necessary in order to achieve a substantial overlap being 
accompanied by charge-carrier derealization. Upon low- 
ering the temperature, obviously the number density of 
MP increases and smaller fields are needed to reach the 
percolation threshold. Also, the larger magnitude of R^ 
indicates a higher degree of inhomogeneity due to the 
larger number of polaronic objects. Since the normal- 
ized measurements collapse onto one curve |3S], a single 
function, see inset of Fig.[3]ja), describes the density of 
magnetic polarons at different temperatures and mag- 
netic fields. 
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